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Abstract 



The Leibniz's rule for fractional Riemann-Liouville derivative is 
studied in algebra of functions defined by Laplace convolution. This 
algebra and the derived Leibniz's rule is used in construction of ex- 
plicit form of stationary-conserved currents for linear fractional dif- 
ferential equations. The examples of fractional diffusion in 1+1 and 
the fractional diffusion in d+1 dimensions are discussed in detail. The 
results are generalized to the mixed fractional-differential and mixed 
sequential fractional-differential systems for which the stationarity- 
conservation laws are obtained. The derived currents are used in con- 
struction of stationary nonlocal charges. 
PACS : 11.30-j 
MSC: 26A33 
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1 Introduction 



In the paper we shall study the properties of the fractional-differential 
equations together with the mixed models containing both fractional 
and standard classical derivatives. 

The fractional analysis describing the fractional integrals and deriva- 
tives is covered extensively in literature (see for example the mono- 
graphies ffl, pi pi 01 an references given therein). Recently these oper- 
ators have found application in various areas of physics. Let us start 
with fractional mechanics describing the nonconservative systems de- 
veloped by Riewe [||, [| who also shows the possible connection be- 
tween the fractional formalism and a problem of classical frictional 
force proportional to velocity. 

The fractional operators emerge also as the infinitesimal generators of 
coarse grained macroscopic time evolutions 
mine fractional diffusion processes 



12, 12, 



|, [lC|, 11] and deter- 



E MM- 

The phenomenological approach to derivation of the stress-strain rela- 
tionships which tends to proper description of the rheological proper- 
ties of wide classes of materials leads to rheological constitutive equa- 
tions with fractional derivatives [17|. 

Next domain is the path-integral formulation of classical boundary 
problems with fractal boundaries used in polymer science. These mod- 
els can be rewritten in form of fractional differential equations. The 
order of the fractional operator is given by the geometry of the bound- 
ary, the space in which the boundaries are embedded and the type of 
random walk process [jD], |l8|| . 

One should also mention the description of wandering processes given 
by the fractional Fokker-Planck-Kolmogorov equation in the fractal 



space-time [19, 2C, 21, 22 1, the fractional generalization of Klein- 
Kramers equation which yields the fractional Raleigh and Fokker- 
Planck models [23, 24, and the fractional equation describing the 



end-to-end distribution of stable random walk where the fractional 
power of the standard Laplace operator is used [18|. 
Finally the fractional operators appear also in field theory where re- 
cently the roots of the wave operator were investigated by Zavada 
[26, In his paper he shows that the Dirac operator is the only 
one from them which can be realized using the standard derivatives. 
When the root of order different from ^ is considered we obtain the 
fractional differential equation. 

Most of the above examples are linear equations with constant co- 
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efficients of mixed type - containing both fractional and standard 
derivatives. As is well-known in classical field theory the conservation 
laws for linear differential systems can be derived using Takahashi- 
Umezawa method [Eql. This procedure has been extended to discrete 



and noncommutative models fl29L 3C , Bl 



Our aim is to show that similar procedure can be applied to frac- 
tional equations in the convolution algebra of functions in order to 
obtain the stationarity-conservation laws which are analogs of conser- 
vation equations known for models from classical differential calculus. 
The explicitly derived stationary-conserved currents are nonlocal ex- 
pressions with respect to this part of space for which the fractional 
derivatives appear in the initial equation. This phenomenon is con- 
nected with the nonlocality of fractional operators as well as with the 
convolution algebra of functions which we introduce to simplify the 
Leibniz's rule in the fractional differential calculus. 
Some of the derived nonlocal currents yield the stationary charges 
which in turn can be converted into nonlocal conserved charges. In 
the present paper we discuss this procedure on some examples and 
then for general case of mixed fractional and differential equations. 
These nonlocal integro-differential equations obey new type of conser- 
vation law which we call stationarity-conservation law. 
In section 2 we review briefly the properties of Riemann-Liouville frac- 
tional integrals and derivatives and show that the Leibniz's rule is sim- 
plified in the algebra of convolution. The new Leibniz's rule produces 
strict requirements concerning the behaviour of the functions in the 
neighbourhood of 0. The next section contains the detailed discussion 
of the derivation the stationarity-conservation law for two examples 
of fractional diffusion: in 1+1 and d+1 dimensions. It is explicitly 
proven that the asymptotic properties of the solutions for diffusion 
equation in 1+1 dimension allow construction of stationary currents 
and stationary charges. Then the currents and charges are converted 
via convolution to conserved currents and charges which are station- 
ary in a strict sense - that means true constant functions. 
Final section includes the general fractional-differential model as well 
as the sequential fractional-differential one. We show explicit con- 
struction of stationary currents, the derivation of the stationarity- 
conservation laws and close the section with discussion of the possible 
stationary and conserved charges. 
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2 Properties of fractional integrals and 
derivatives 



2.1 Riemann-Liouville fractional integral 

Let us recall the definition of Riemann-Liouville fractional integral 
|I[ 0, U used widely in the literature dealing with fractional calculus: 

Definition 2.1 Let Rev > and let f be piecewise continuous on 
(0, +oo) and integrable on any finite subinterval of [0, +oo). Then for 
t > 

ZV7(i):=^7X !\t-sy- l f(s)ds (1) 



o 



is £/ie Riemann-Liouville fractional integral of f of order v 



We notice that the above definition includes the operation of Laplace 
convolution, namely it can be written as: 

D^f(t) = * f{t) = f * *_„(i) (2) 

where we have denoted <fr- u (t) = f^y^ -1 - 

Now we are interested in the properties of the fractional integral con- 
nected with the composition of the integrals with respect to the same 
coordinate. The answer is the generalization of the Dirichlet's integral 
formula for continuous function which is called in fractional calculus 
the composition rule ^, [J : 

D- u D-^f(t) = D-^ + ^f{t) = D-^ L D- U f{t) (3) 

for Refi, Rev > and for any function / piecewise continuous on 
[0,+oo). 

Let us now present the known forms of Leibniz's rule for integral (|l|): 



— V 



D-»(f.g) = '£(j)D-»-3f-Dlg (4) 

where / and g are real analytic functions on [0, +oo). This rule was 
generalized by Osier [Q, [3^, |33|, 34] who obtained the following forms 
of Leibniz's rule: 

+00 -p/ -| \ 

(5) 
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D~ u (f-a)= ( + °° r ( v + l ) D -u--y-\ f.n7+A j\ 

[I9) r(- I /- 7 -A + i)r(7 + A + i) jU ; ^ 5dA 

(6) 

where 7 is an arbitrary complex number. 

We shall not discuss the convergence of the series in (|5|) and of the 
improper integral in (||). Let us notice however that when the algebra 
of functions is defined by standard point-wise multiplication as in the 
above formulas all versions of Leibniz's rule are very complicated. 
Thus we propose to investigate the algebra of functions with multipli- 
cation defined via Laplace convolution: 

f*g(t):=f t f(t-s)g(s)ds (7) 
J 

As is well known this multiplication is associative and commutative . 
The neutral element is the Dirac ^-function. Let us prove the following 
Leibniz's rule for fractional integral (|l]) and multiplication defined by 

D-"(f*g) = (D-Wf)*D-Vg (8) 

where Rev > and 7 a complex number fulfilling inequality Re(v — 
7) >0. 

The new Leibniz's rule is implied by the composition rule (E3) and 
properties of convolution which defines the fractional integral @ and 
algebra of functions (|?]): 

D-"(f * g) = D-7-(-7)(/ * g ) = D^D-^\f * g) = (9) 

(/ * g) * 3>_ ( „_ 7) ) * $- 7 ) = (/ * $_(„_ 7 )) *(g* $- 7 ) = 

= D-^-t) f * D~~i g 
where i?ez^ > and i?e(z/ — 7) > 0. 

The derived formula (^) is similar to the multiplicity properties of 
the transformation operators in the discrete [29] and noncommutative 
[30, [H]] differential multidimensional calculi for standard product of 
functions: 

q(f-g) = (tif)-($g) (10) 

The multiplicity property of the fractional calculus (||) leads to the 
following redefinition of the integral of order v. 

DTVW := (Dt U ~ !)/(«) = / * (*-„ " W) (11) 



6 



The new operator D obeys the following Leibniz's rule in the algebra 
defined by convolution product (0): 

T>r (/ * 9) = (C 7 /) * + (Dp7) * g (12) 

or its symmetric form: 

Vt»U *<?) = /* A~ 7 <? + iv;^f) * Cg (13) 

where for given 1/ the 7 is a complex number fulfilling conditions: 
Re~/ > 0, Re(v — 7) > 0. As we have noticed the analogy between 
the action of the fractional operator D in the algebra of convolution 
product (^) and the transformation operator £ in the dicrete and non- 
commutative algebra (|l0|) we shall use in the sequel the notation £ for 
the "old" fractional integral ([!]): 

C a = (14) 

The above Leibniz's rules are implied by the properties of the con- 
volution and the composition rule (||). If v and 7 fulfill the above 
restrictions we obtain: 

vr(f*g)= (15) 

/ * g * ($_„ - 6) = f * g * ($_ 7 _ (i ,_ 7) ± $_ 7 - 5) = 
/ * g * $_ 7 * (*_(„_ 7 ) - 6) + / * g * ($_ 7 - 5) = 
(/ * $_ 7 ) *g* ($_ (v _ 7 ) - 5) + / * ($_ 7 - 5) * g = 

{C 1 f)*V^ v -i ) g + {Vpf)*g 
The proof of the symmetric form of the Leibniz's rule (13) is analogous. 



2.2 Riemann-Liouville fractional derivative 

The operator known as the Riemann-Liouville fractional derivative 
HI H|, m is defined using the fractional integral (jl]) : 

Definition 2.2 Let m < Rev < m + 1, t > 0. TTie operator given by 
formula: 

A" := I ^ ) A" (m ^ +1) /W (16) 



/or functions for which the improper integral on the right-hand side of 
( |7^ j is convergent is called the Riemann-Liouville fractional derivative 
of order v . 
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Let us notice that the functions from the domain of the operator 

form the subset in the set of functions from definition 2.1. It is well- 
known fact that this class consists of finite sums of functions of the 
type: 

oo 

t A E«^ ( 17 ) 

or 

oo 

t x ln(t) a k t k (18) 

k=0 

where ReX > — 1 and the series have a positive radius of convergence. 
Contrary to the fractional integrals the derivative ( |l~6l) cannot be ex- 
pressed using only convolution. The formula includes the classical 
derivative and looks as follows: 



d 



m+l 



D»f{t) :=[-) (f**u-m(t)) (19) 



with the function § v -rr, ~ 



r(m+l-f) " 

We expect the fractional derivative to obey the composition rule anal- 
ogous to the one for fractional integral. In fact [|l], 0] the following 
formula which generalizes (J3j) is valid: 

DZD?f = Dp»f (20) 

provided: 

• v arbitrary, \i < X + 1 and the function / is of the type described 



by (13,18) 



v arbitrary, fi > A+l and = k = 0, ...m—1 for the function 



/ of type (17,18) where m is the smallest integer greater or equal 
to Re/i. 

The above formula shows that the fractional derivatives of different 
orders do not always commute as it is the case with the fractional 
integrals. 

The Leibniz's rule for fractional derivative has the form {Rev < n— 1) 
§0||||g]: 



D tf ■ 9(t) =£(*) 9 {k) ■ Dr h f(t) ~ R u n (t) (21) 
k=0 v 1 



<s 



when the function / is continuous in the interval [0, t] while g has n+ 1 

continuous derivatives in [0, t]. 

The remainder R n is the integral expression: 

Kit) = -^i— _ j\t- S y- l f(s)dsf s g( n+1 Hu)(s-urcL; (22) 

If the above remainder goes to for n — ► oo the Leibniz's rule (21) 
can be written for analytic functions in the form of series: 

00 / \ 

Dtf-9 = Y, [ k )f {k) -D v t - k 9 (23) 
k=o ^ ' 

Again the form of Leibniz's rules for the algebra defined by point-wise 
multiplication of functions is complicated. 

We propose to use the algebra of convolution ([?]). The following stat- 
ment is valid for the new algebra of functions: 

Lemma 2.1 Let m < Rev < m + 1 and the function g be piecewise 
continuous in (0, +oo). If the function f is a finite sum of functions 
of the type ( [7^ , and fulfills the condition: 

lim /( fc )*^_ m = 
for k = 0, l..,m then the following rule holds: 

DUf*g) = (DU)*g (24) 

Proof: 

We use the well-known rule for differentiation of an integral depending 
on a parameter with the upper limit depending on the same parameter: 

^ [ t F(t,s)ds= [ t ^Ads+ lim F(t,s) (25) 
dt Jo Jo at s^t-o 

and from it follows for < Rev < 1: 

D v t U*9) = J t (f*9*®») = J t (f*®»*9) = (j t (f * *9 = {D v t f)*9 

(26) 

provided: 

lim / * $Jt) = (27) 
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Thus when the assumptions are fulfilled the formula ( |24| ) is valid. 
For m < Rev < m + 1 we apply the rule (p5|) m + 1 times: 



j v m+1 



^ x m+1 

— J [(/ * ^_ m ) * 5 ] 



(/**i/-m)] *5 



m+1 



(/ * 



*5 



and arrive at the conditions: 



lim t ^o+o / * $u-m{t) = 
lim^o+o /' * ®v-m(t) = 

lim^ +o/ (,n) *^-m(t) = 



(28) 
(29) 
(30) 
(31) 



which are fulfilled by assumption. 

The above set of right-sided limits determines the behaviour of the 
function / in the neighbourhood of t = 0, namely f(t) ~ t™ with (3 a 
complex number fulfilling the condition: Re/3 > — 1 + Rev . 



The symmetric version of the formula (|24| ) follows from the commu- 
tativity of the Laplace convolution. 

Corollary 2.2 Let m < Rev < m + 1 and functions f and g are 
piecewise continuous in (0, +oo). // both functions f, g are finite sums 
of functions of the type §T\ , and both of them obey the assumptions 
from Lemma 2.1 then the following rule holds: 



DUf * 9) = &{D»f) *g + Q.-/3)f* (D»g) 



(32) 



for (5 G [0, 1] 



The above lemma together with the composition rule ( |20[) yields the 
analog of the property (||) for Riemann-Liouville fractional derivative: 



Corollary 2.3 Let Rev > and the function f * g obey for certain 
7, fulfilling Rej > and Re(v — 7) > 0, the assumptions of the 
composition rule f^j. If function f fulfills the conditions from Lemma 
2.1 for 7 and the function g the corresponding conditions for v — 7 
then the following formula holds: 

DUf * g) = Dt~ 7 Dj(f * g) = (Djf) * D^g (33) 
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Analogously to ( ]TT| ) we can introduce the new differintegral operator: 

2#/(t) := W - I) fit) (34) 

The Leibniz's rule for the introduced differintegrable operator of pos- 
itive order v is similar to the one known from the discrete and non- 
commutative calculus [|^, |3(| j31~l : 

W * 5) = (A 7 /) * 9 + (C 7 /) * Vt^g (35) 
* 5 ) = [Vjf) * C~ J 9 + / * T>^g (36) 

where we use the notation: 

C = DJ (37) 

and v, 7 together with functions /, g fulfill the conditions from Lemma 
2.1. 

2.3 Riemann-Liouville partial fractional deriva- 
tives 

Let us extend the formalism introduced in previuos sections to multi- 
dimensional case. We shall study the stationarity-conservation equa- 
tions for some fractional partial differential equations and derive for 
them the explicit form of stationary currents connected with their 
symmetries. We assume that in the equation both types of derivatives 
can appear - fractional with respect to to a subset of coordinates and 
classical - continuous ones with respect to the rest of coordinates. 
Thus the question arises how to define the multiplication of func- 
tions. We propose to use the multidimensional Laplace convolution 
when initial equation contains only Riemann-Liouville fractional par- 
tial derivatives of the form: 

D a k k f(x) := (38) 

w -A T (d* k ) mk+1 f V " s)- a ^f(x + (s - x k )e k )ds 

T{m k + l-a k ) Jo 

where m k < Rea k < m k + 1. The upper index in the formula denotes 
the fractional order of the partial derivative while the lower one says 
that it was taken with respect to coordinate x k . 

Let subset of coordinates in our n-dimensional model 

for which the fractional partial derivatives (pS|) appear in the equation. 
Then we define multiplication of functions as follows: 
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Definition 2.3 The algebra of functions is defined by the multiplica- 
tion formula: 

f * 9(x) := (39) 
y o -y o f -^sieij g + 22(si - xi)eij ds 1 ...ds m 
where {e{] k = S lk . 

Similarly to the one-dimensional case the multiplication (|39|) is asso- 
ciative and commutative. 

In the above algebra of functions the Leibniz's rule ( |32| ) given by 
Corollary 2.2 is valid for functions fulfilling the respective assump- 
tions concerning their behaviour at x k = 0: 

D a k k f * g = f3 k (D^f) *g+(l- k )f * D?g (40) 

with (5 k £ [0, 1] for k = 1, ...,m. 

For classical derivatives acting by assumption in directions j = m + 
1, ...,n we obtain for convolution (^) the standard form of the Leib- 
niz's rule: 

d j (f*g) = (d j f)*g + f*d J g (41) 

Similarly to the one-dimensional case investigated in the previous sec- 
tion we can introduce also the partial differintegral operators of posi- 
tive order for functions fulfilling suitable conditions: 

V^f(x):=D^f(x)-f(x) (42) 

These operators obey the Leibniz's rule for functions multiplied ac- 
cording to (|3T 



K k (f * 9) = {Vff) * g + (C7) * K k ~ lk 9 (43) 
V a k k (f*9) = (Vl k f)*Qt lk 9 + f*V a k ^g (44) 
d j (f*g) = (d j f)*g + f*d j g (45) 

where k = 1, ..,m and j = m + 1, ...,n and the function / obeys the 
conditions of Lemma 2.1 for the fractional order of the derivative 7& 
while the second function g respectively fulfills these conditions for 
a k - 7fc- 

The first two formulas are the symmetric forms of the Leibniz's rule 
for fractional derivatives and the last one is standard Leibniz's rule 
for partial derivatives but taken in algebra of functions defined by 
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multiplication (p9|). 

Now we can apply the properties of multiplication (|3^) and fractional 
differentation in construction of the stationarity-conservation laws and 
conserved charges for some partial fractional equations. 



3 Examples 

3.1 Fractional diffusion equation in 1+1 

Let us recall the fractional diffusion equation discussed in || H, 14 1 : 

D?4>(x, t) = X 2 d 2 ^(x, t)+ct>(x,0) \ - (46) 

l (1 — a) 

where t > 0, x £ R and < a < 1 describes the process of ultraslow 
diffusion while the value 1 < a < 2 is used for intermediate processes 

111,0- 

Let us focus on the case of ultraslow diffusion. The operator of the 
equation contains both types of derivatives: fractional with respect to 
time and standard for the spatial dimension: 

A(D?,d x ) = D?-\ 2 d 2 x (47) 



The product of functions for this model is defined according to ( |39[) 
and looks as follows: 

f*g{x,t)= [ f{x,t-s)g(x,s)ds (48) 
Jo 

Using the properties of the new multiplication and of the fractional 
derivative (32) we construct the operator T with components: 

= X 2 d x - \ 2 d x T t = 2 (49) 

Then the current: 

j x = <j)'T x * <\> = ft\ 2 d x *<!>-<//* \ 2 d x 4> (so) 

J t = <j)'T t * (f> = 2(f)' * <f> (51) 

obeys the stationarity-conservation equation for t > in the area 
cj)(x,0) = <jf{x,Q) = 

8 X J X + D?J t = (52) 
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provided the function is the solution of initial equation ( f46| ) while 
<j)' solves its conjugation: 



A(-£ff , -d x )cf>'(x, t) + ef>'(x, 0) \ a = (53) 

1 (1 — a) 



Before passing to the proof of the stationarity-conservation law ( p2[ ) 
we shall discuss the existence of solutions of diffusion equation and of 
its conjugated form with required properties around t = 0. 
Let us recall the form of general solution of the equation (^) (lj] : 



<f>(x,t) = / dyG a (x,y,t)<t>(y,0) (54) 

J — oo 

where G a is the fractional Green's function of the following form: 

/*oo /*oo 

G a (x,y,t) =t~ a dzE a (t- a z)G(x,y,z) = dvE a (v)G(x,y,t a v) 
Jo J 

(55) 

with the function G(x,y,z) = G(\ x — y \,z) being the standard 
Green's function: 

1 \x-y\ 2 

G(\ x-y \,z) = -j=e ^ 

V 47TZ 

and E a denoting the Mittag-Leffler function [Q, Q . 
Taking into account the asymptotic properties of the function G we 
conclude that the solution <j) behaves in the neighbourhood of t = as 
the power function t~? . The solution cf>' of the conjugated equation 
has a similar form so its behaviour for t — > is the same as of the 
considered solution <p. 

This fact implies that at least for < a < | the assumptions of 
the Lemma 2.1 are fulfilled therefore we can use in the proof of the 
stationarity-conservation equation (52) the Leibniz's rule for fractional 
derivative Df given in formula (|3^), 

Let us check the conservation law explicitly applying the Leibniz's rule 
OH) with = \: 

8 X J X + DfJ t = (56) 
d x (V A 2 ^ * - <P' * X 2 8 x <p) + Df (2<p' * <p) = 
\ 2 {d 2 J) *4>-<p'* X 2 d 2 x <j ) - (-AV) * <f> + 4>' * D?<t> = 
- '(-D? - X 2 d 2 x ) 0'] * + 0' * (A° - A 2 d 2 ) = 
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We have omitted the terms depending on initial values 4>{x, 0) and 
4>'(x,0) as we expect the rule (|52|) to be fulfilled in the area where 
c/)(x,0) = <f>'{x,0) = 0. 

Having obtained the general form of stationary current ( |50| , 51 ) we can 
discuss the possible symmetries of equation (|46|) which can be used 



in construction of different solutions of the initial diffusion problem. 
The set includes spatial momentum P x = d x as this operator commuts 
with the operator of diffusion equation fl46|). 

The stationarity-conservation laws for currents including new solu- 
tions are fulfilled for transformed solution P x (f> in the area d x (j)(x, 0) = 
4>'(x,0) = 0. The simplest possible choice of initial value for solutions 
of diffusion equation and of its conjugation is = - = S(x) 

with cpo and cp' arbitary constants. 



In this way we arrive at the stationary (for x ^ in the sense of (52)) 
current connected with symmetry of the fractional diffusion equation: 



6T T * P T 



TX 



*T t * P x <, 



(57) 



The stationarity-conservation equation (52) can be reformulated us- 
ing the definition of Riemann-Liouville derivative so as to obtain the 
standard conservation equation namely: 



d x j' x + d t 4 = o 



(58) 



which is fulfilled for i / and the components of the new current 
look as follows: 



V - J 



J't = Jt*$a 



1 



r 1 - a 



-Jt * t- a 



(59) 



Following the classical field theory the time-components of the con- 
served currents J and J' yield the charges: 



Q = IZodx Jt 
Q' = JZodx J[ 



(60) 
(61) 



The respective derivatives of the above charges are determined by the 
boundary terms for time-components of the currents and the initial 
conditions for solutions <b and d>': 



D?Q= lim \\ 2 {d x 4>')*<l>- AV *d x (j) 



+ 



(62) 



— lim 

x — >~oo 



A 2 (9^')*0- \ 2 <P'*d x ct> 



+ 



15 



T(l-a) rv ' ' r v ' ' ru r(l-a) 

|q' = (63) 



From he general form of solutions ( pl[ ) we obtain for initial condition 

4>(x,t = 0) = <po5(x): 

0(0, t) =0 o G Q (O,O,t) 

and for conjugated equation: 

0'(O,t) = -0' o G a (O,O,i) 

Due to this property of the solutions and the commutativity of the con- 
volution the last terms in the above formulas cancel. The first parts 
vanish by the asymptotic properties of the Green's function which de- 
cays exponentially together with its spatial derivative for large x. 
Thus the explicit expressions for charges (|60,|61) produce the station- 
ary function Q and constant function Q' connected with the station- 
arity law and conservation law of the diffusion equation in 1+1 dimen- 
sions: 

DfQ = j t Q' = (64) 

3.2 Generalized fractional diffusion 

Let us now extend the dimension of the space-like coordinates to d. 
We shall consider the equation known as the generalized fractional 
diffusion problem [jrj, 35]: 

D?4>(x, t) = CA<t>(x, t) + <f>(£, 0) * - (65) 

L(l — a) 

with t > 0, < a < 1 and A the Laplace operator in d-dimensional 
Euclidean space. 

We again check the properties in the neighbourhood of t = of the 
solutions. 

The solution for arbitrary initial condition generalizes the formula ( |54| ) 
used in the previous section: 

POD 

<f>(x, t) = / d d yG a (x, y, t)0(y , 0) (66) 
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where G a is the fractional Green's function of the following form: 

roc 

G a {x,y,t) = r a dzE a (t- a z)G(x,y,z) (67) 



with the function G(x,y,z) = G(\ x — y \,z) being the standard 
Green's function: 

G(\ x — y\,z) = (4irz)~ 2 e 4T - 
and E a denoting the Mittag-Leffler function. 

Taking into account the fact that this function for < a < 1 is an 
entire function and vanishes exponentially for large positive values of 
argument we conclude that the solution cj> behaves in the neighbour- 
hood of t = as the power function t~ a . Similar argument applies to 
the solution of the conjugated equation (j)' given below (pi). 
The product of functions given by ( |39|) has in d + 1-dimensional case 
the following explicit form: 

f*g(x,t)=f f(x,t-s)g{x,s)ds (68) 





The number of the components of the operator T and of the current 
J is now d + 1 while the form of the space-like and time-like parts is 
identical to the ones obtained for the modified Nigmatullin's diffusion 
equation. 

The operator T given by: 

Ti = Cd Xi - C8 Xi F t = 2 (69) 

can be applied in the construction of the current: 

Ji = cP'Ti *<j) = 0>Cd Xl *(b-cp'* Cd Xi 4> (70) 
J t = (f,'r t * <f> = 2(j)' * (71) 



where (ft solves the initial generalized diffusion equation ( pq ) and (j)' 
its conjugation : 

A(-D?, -d Xl , -d Xd )<j>'{x, t) + <t>\x, 0) * - = (72) 

1 (1 — a) 

The current (70, 7^) obeys the stationarity-conservation equation: 

d 

J2d Xi Ji + D?J t = (73) 

i=l 
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for x ^ provided the solution of diffusion equation ( pq ) with the 
initial condition (j)(x,0) = 4>o5(x) is taken and for the conjugated 
equation the initial condition (f>'(x, 0) = 4>'q5{x) is considered. 
The proof of the above conservation law is analogous to the one pre- 
sented in the previous section for the 1 + 1 diffusion equation. The 
essential feature in the proof are the asymptotic properties of the solu- 
tions (f) and 4>' in the neighbourhood of t = 0. Similarly to the previous 
case they allow us to apply the Leibniz's rule (32) for x ^ at least 
when < a < ^. 

The set of symmetry operators for d + 1 case is much wider as it 
contains not only momenta: 

Pi = d Xi (74) 

but also the angular momentum with respect to the space-like dimen- 
sions: 

= Xid X] - Xjd Xi (75) 

where i,j = 1, d. 

As the symmetry operators transform solutions of ( |65| ) into solutions 
with the same properties around t = we can use them in construction 
of the stationary-conserved currents: 

Jf = 4>'Ti * 5(f) = if/Cdxt * &<t> ~ 4>' * Cd Xi 5cP (76) 
Jf = cj)'Tt *5<j> = 20' * 5<j> (77) 

where 5 denotes one of the above symmetry operators of the equation 

dH). 

The currents ( [7^j77| ) can be transformed into the components J' sim- 
ilarly as in the case of the 1 + 1 fractional diffusion. Taking the com- 
ponents: 

J'i=4 4 6 = J$*$ a = 1 , J?*t~ a (78) 

1 (1 — a) 

we obtain the conservation law for d + 1 fractional diffusion process: 

</ 



J2d Xi Jf + dtJ't S = (79) 
i=i 



valid for x ^ 0. 

Finally the derived set of stationary and conserved currents yields two 
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sets of charges indexed by the symmetry operators 5 € {Pi, Mij} i, j = 
l,...,d: 



which are respectively stationary and conserved functions of time: 



4 Stationarity - conservation laws for 
some fractional partial equations 



In previous sections we have obtained the stationarity law and con- 
servation law for some examples of partial fractional differential equa- 
tions. 

In the multidimensional case of diffusion equation the general solu- 
tion allows the explicit construction of the current which obeys the 
stationarity-conservation law in the area of space where the initial 
conditions vanish both for solution of diffusion equation and for its 
conjugation. 

We have shown for that the stationary currents yield stationary charges 
which can be converted to the conserved ones. 

The discussed example shows that the construction of possible charges 
stationary or conserved is connected with the asymptotic properties 
of solutions around and for | x |— > oo. 

In the sequel we shall discuss the general construction of stationarity- 
conservation laws assuming that regular in the sense of Lemma 2.1 so- 
lutions of the respective fractional differential equations exist at least 
in certain area of space. 

4.1 Mixed fractional differential and differen- 
tial partial equations 

Let us consider now the general equation which contains the fractional 
and differential parts of the following form: 



Q S = C d d x Jf 
Q' S = f™d d x J* 



(80) 
(81) 



DfQ 5 = 
iQ' s = 



(82) 
(83) 



A{D,d)(f>= [A(D] + A(d)](f> = 



(84) 
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+A U = 



The introduced equation has constant coefficients (we admit also con- 
stant matrices) and generalizes the studied examples of fractional dif- 
fusion. We shall study the construction for the homogenous form of 
the equation remembering that the addition of the initial terms similar 
to the ones discussed previously restricts only the area of application 
of the stationarity equation and does not change the general con- 
struction. We assume that for given variables x%, ■■■,x m the equation 
includes only fractional derivatives in A(D) while for the remaining 
coordinates x m+ i, x n only partial derivatives appear in the opera- 
tor A(d). 

To derive the stationarity-conservation law we shall use the Takahashi- 
Umezawa method [^8| for the differential part A(d) and the fractional 
Leibniz's rule ( |32|) for the part A(D) containing fractional operators. 
As we know from the discussed examples each direction of the space 
yields the component of the current which for coordinates X\, .-,Xk is 
given by the T operator of the form: 

f k = 2A k (85) 

while for the part j = m + 1, n we obtain [|28|: 



r,= EEV,,(-9 )-(-d (86) 

1=1 k=l 

It is the well-known fact that for an arbitrary pair of functions / and 
g the operator T fulfills the equality: 

n _ 

]T ^/ * T jg = -/A(- d)*g + f* A(d)g (87) 

j=m+l 

where the multiplication is given by the convolution (^) and A(— d) 
is the conjugated operator for A(d) acting on the left-hand side. 
The above property of the T operator together with the Leibniz's rule 
(|32l ) for fractional derivatives (taken with parameters (3^ = \ k = 
1, ...,m) implies the following proposition to be valid: 

Proposition 4.1 Let the function <p be an arbitrary solution of the 



equation (84) and let (j)' solve the conjugated equation: 
<£'A(- D, - d) = <t>'[A{- D) + A(- d)\ = 
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-E A ^ +EV...m ; (-9 ) + a =o 

fe=i i=i / 

TTien the current given by the components: 

J k = $ * A k (fi + (j)'A k * (j) k = l,...,m (89) 
Jj = (j)' * j = m + l, n (90) 

fulfills the stationarity- conservation equation: 

m n 

J2 D k k Jk+ E ^ J i = ° (si) 

fe=l j=m+l 

provided the solutions (j) and <p' fulfill the conditions of Lemma 2. 1 in 
the neighbourhood of x k = k = 1, ...,m. 

Proof: 

We check the law ( plf ) explicitly: 

m n 

E^ + E ^ = 

k=l j=m+l 
m n 

E^ fc (V '*h<t> + <P % *</>) + E 9 , '(0 / *r i 0) = 

fc=l j=m+l 
m m 

E(^r^')Afc * + E>' * A*l£ h - 0'A(- a) * + 0' * A(0)0 = 
fc=l fc=l 

-^'A(- D, - d) * ^ + (j)' * A(D, d)4> = 



Thus for every equation of the form (|84j ) we can produce exact form 
of the stationary-conserved current provided the initial equation and 
its conjugation have solutions which fulfill the asymptotic conditions 
at Xfc = k = 1, ...,m which allow the application of the Leibniz's 
rule for fractional partial derivatives Du k . 

The stationarity-conservation equation (91) can be rewritten in the 
form of the standard conservation law for modified components of the 
above current (m k < a k < rn k + 1 k = 1, m): 

J k = (d k ) m «(J k * k <S>a k -m k ) k = l,...,m (92) 
Jj = Jj j = m+ 1, n (93) 
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where the convolution * k is given by the formula: 

/ *k g(x) = (94) 

r x k 

\ f(x- se k )g(x + (s - x k )e k )ds k 
Jo 

The new current J' obeys the conservation law: 

n 

£<9<J/ = (95) 
1=1 

4.2 Mixed fractional sequential and differen- 
tial partial equations 

In the previous construction we have considered the fractional part 
of the operator including only the first power of the corresponding 
partial fractional derivatives while in the differential part we have 
taken an arbitrary polynomial of partial derivatives. Let us extend 
the derivation of the stationarity-conservation laws to the general case 
containing both polynomial of fractional derivatives and polynomial 
of classical partial derivatives: 

A(D,d)<f> = [A(£>) + A(9)]0= (96) 




The derivatives with respect to the coordinates x±, x m are the frac- 
tional Dp* where the upper index denotes the fractional order and the 
lower one the respective partial direction. The part depending on frac- 
tional derivatives has now the form of partial sequential fractional op- 
erator generalizing the sequential operator for one-dimensional space 
H). The coefficients A and A are again constant matrices or numbers. 
As the derivatives with respect to different coordinates do commute 
both types of coefficients are fully symmetric with respect to the per- 
mutation of the set of indices. 

To obtain the T operator fulfilling the equation (j87|) we again use the 
Takahashi-Umezawa method for the differential part A(<9) and get the 
components Tj as given by ( p6| ) whereas for T we have: 

M-i j 

T k = 2 £ J2A kpi ___ Pj (- Dpi )...(-D pi )D^...D^ (97) 

j=l 1=1 
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It is easy to check the analog of the formula ( |87| ) for the operator f : 

m 

£ (/ * f k g) = -/A(- 5) * g + / * A(D) fi (98) 

k=l 

for an arbitrary pair of functions / and g allowing the use of the Leib- 
niz's rule (^) together with their fractional derivatives Dp}** ...D p !j g 

andf(-D pi )...(-D pi ). 

All the above calculations yield as a result the following proposition 
which describes the explicit construction of the stationarity - conser- 
vation law for linear sequential fractional-differential equation (|96|): 



Proposition 4.1 Let the function <f> be an arbitrary solution of the 
equation ( 96 ) and let <j>' be a solution of the conjugated equation in the 
form: 

= (j)'A(- D,-d)= (99) 

- «-ai <-a k * _ W \ 

qEV., i (-D w )...(-J}J|^A,.., 1 (-a )...(-d )+A 
\fc=l 1=1 ) 

Then the current with the following components: 

J k = <//*t k </> k = l,...,m (100) 
Jj = 4>' *Tj(t> j = m + l,...,n (101) 

obeys the stationarity- conservation equation: 

m n 

J2 D k kJ k+ E ^ J i = ° ( 102 ) 

fc=l j=m+l 

provided the solutions (j) , <p' together with their derivatives appearing 
in the formulas for components ( 10C\) fulfill the conditions of Lemma 
2.1 in the neighbourhood of x k = k = 1, ..,m 

Proof: 

We use the properties of the solutions and of the operators T and T 
and obtain: 

E]= m+ i d j Jj = E]= m+ i d J (</>' * T j(p ) = -<f/A(- d)*<f> + 4>'* A(d)<t> 
ET=i D a k k Jk = E2Li U' * F k A = -<j>'~A{- D)*<t> + <t>'* ~A(D)cp 
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Thus the left-hand side of the stationarity-conservation formula is of 
the form: 

m n 
k=l j=m+l 

-0' (A(- 5) + A(- d) + A ) * 4> + <(>' * (A(D) + A(d) + A ) = 
and vanishes on shell. 



We can rewrite the stationarity-conservation law to have the conser- 
vation law connected with the equation (^). To this aim we apply 
the definition of the Riemann-Liouville fractional derivative (^) . The 
modified components of the current have the form similar to the one 
derived in the previous section (m& < a k < + 1 k = 1, m): 

J' k = (d k ) mk (Jk*k^a k -m k ) k = l,...,m (103) 
Jj = Jj j = m + l,...,n (104) 



with the convolution *k given by (j& 
They obey the conservation law: 



E^' = ° ( 105 ) 



i=i 

4.3 Stationary and conserved charges for mixed 
fractional-differential models 

Following the results obtained for fractional diffusion we shall apply 
the derived stationarity-conservation law in construction of stationary 
charges. 

Two cases should be considered: when the time-derivative in the op- 
erator of the equation is a fractional and when it is standard partial 
one. 



Let us assume that the time-derivative in equations (84,96) is a frac- 
tional one. Integrating the time-component of the current fulfilling the 
stationarity-conservation equation ( pl|jl02 ) we arrive at the charge: 



dx J t (x,t) (106) 

I Ft™- 1 

which is a stationary function of order at which also determines the 
order of the fractional time-derivative: 

D^Q = (107) 
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provided the respective boundary terms vanish. For components Jj j = 
m+ 1, n it means that they vanish at the infinity in the given direc- 
tions while for components J k k = 2, ..,m the asymptotic condition 
has the form: 

lim (&) m >>(J k * k *a h -m h ) = (108) 
\x k \-*oo 

where m k < a k < m k + 1. 

The second possibility is the model with standard time-derivative. 
Then the charge: 

Q= I dx J t (x,t) (109) 

is a strictly stationary function of time that means it is a true constant 
function: 

d l Q = (110) 

when the asymptotic conditions for respective components of the cur- 
rents are fulfilled: 

lim Jj=0 j =m + 2,...,n (111) 
lim (d k ) mk (J k * k <S> ak - mk ) = k = l,...,m (112) 

\x k \^oo 

(113) 



The exact form of the symmetry algebra of the equations (84,96[) vary 
for different examples. Let us however notice that it includes for all 
of them the momenta: 

Pk = D^ k k = l,...,m (114) 
Pj = &> j = m + l,...,n (115) 

as they commute with the operator of these equations. 
However if we propose to use the above momenta in derivation of con- 
served currents and charges we must additionally assume the regular 
behaviour of the W(D)P k cp and W{D)Pj<p functions in the neighbour- 
hood of with respect to the xi, x m coordinates (W{D) denote the 
polynomials of fractional derivatives appearing in the formula for V 
operator). 

When this assumption is fulfilled the stationary-conserved currents 
look as follows: 

Jl = (j)' *f k 8<t> k = l,...,m (116) 
J 6 . =(//* TjScf) j = m + 1, .., n (117) 
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where the operators T and T are given explicitly in previous sections. 
In this case we have the family of stationary (or respectively conserved 
charges) depending which of the considered two cases apply to our 
model. They have the following explicit form: 



Q s = I dx <j>' *t t 5<j> (118) 



for the case where the time-derivative is fractional and for the standard 
time derivative we have: 

Q 5 = [ dx <j>'*T t 5<f> (119) 
where d is one of the momentum operators given in ( |iT^ , [lT5| ). 



5 Conclusions 

We have discussed the Leibniz's rule for the algebra of Laplace con- 
volution of differintegrable functions. 

The derived procedure for construction of the nonlocal stationary 
currents applies to fractional differential linear equations including 
Riemann-Liouville fractional and classical derivatives provided there 
exist the regular solutions of initial and conjugated equation. Similar 
method is being investigated also for Weyl fractional derivatives for 
algebra of functions defined by Fourier convolution. It seemes that 
it can be extended to models with fractional derivatives defined via 
generalized functions approach as well. 

For the general case we have extracted the explicit form of stationary- 
conserved currents assuming the regularity of solutions in the neigh- 
bourhood of 0. It was shown that the stationary currents are con- 
nected with the conserved ones. Both types of currents produce 
charges: namely stationary currents yield the stationary charges and 
respectively from conserved nonclocal currents we obtain integrals of 
motion. 
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